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We present an experiment of preparing Werner state via spontaneous parametric down-conversion 
and controlled decoherence of photons in this paper. In this experiment two independent BBO 
(beta-barium borate) crystals are used to produce down- conversion light beams, which are mixed 
to prepare Werner state. 
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I. INTRODUCTION 



Entanglement is not only one of the most striking features of quantum mechanics, it also plays a crucial role in the 
field of quantum information Entangled states are important resources for most applications of quantum 

information such as quantum key distribution fl , superdense coding [0 , quantum teleportation || and quantum error 
correction [ffl] etc. Although the best performance of such tasks requires maximally entangled states (Bell states), 
the decoherence effects due to the environment make the pure entangled state into a statistical mixture and degrade 
quantum entanglement in the real word. For a practical purpose, a purification scheme may be applied to the degraded 
entanglement ||J|. 

One of the most important degraded Bell states is Werner state [|l0| , A Werner state in 2 x 2 system takes the 
following form j^j: 

1 — F AF — 1 

Pw = — ^— h + — ^— |*-)(*-|. (1) 
where /„ denotes the n x n identity matrix and \^~) is the singlet state of the four Bell states 

|a>±) =-L(|oo)±|ii)), 
I**) =-^(|01>±|10». 

The Werner state pw is characterized by a single real parameter F called fidelity. This quantity measures the overlap 
of Werner state with a Bell state. In the case where F ^ 1/2, the state is separable, and thus has no entanglement 
to recover or maintain. And a Werner state with F > (2 + 3^) /8 ~ 0.78 violates the Clauser-Horne-Shimony-Holt 
(CHSH) inequality 0,[2|Jl|. 

There are many theoretical investigations on Werner state JTJ], and the Werner state plays an important role 
in entanglement purification [pi, nonlocality Jl(J, entanglement measure |15| and etc. However, there is no report 
of experimental realization of Werner state to date. At present, the most accessible and controllable source of 
entanglement is obtained from the process of spontaneous parametric down-conversion in a nonlinear crystal. In this 
paper, we present an experimental preparation of bi-photon Werner state via spontaneous parametric down-conversion 
and the prepared state is in the form 

(f w = x\*-){*-\ + {\-x) |, (2) 

which can be transformed to pw by a local unitary transformation of a x ® I, where a x is one of Pauli operators and 
/ is the identity operator. 
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There are previous works on preparing mixed sing le photon or bi-photon state like in Ref. (Hg^S A. G. 
White et al. |20| have reported an optical two-qubit source that can produce a wide range of states in the form 

p = x |$+) ($+| + (1 - x) |01) (01| . 

C. Zhang [g2| has proposed a theoretical protocol to produce an arbitrary two-bit mixed state by using beam splitters 
with variable polarization transmission coefficients and single-mode optical fibers. However, in this paper, we propose 
a different way to prepare Werner states conveniently in experiment. 

This paper is organized as follows. In Sec. II, the experimental set-up is described. The experimental results is 
given in Sec. Ill and the discussion and summary is given in Sec. IV. 



II. EXPERIMENTAL SET-UP 



In our experiment, the Werner state is bi-photon's polarization state, which is obtained by mixing a complete mixed 
state and an entangled state. 

The experimental set-up is depicted in Fig. 1. The 1.7-mm-diam pump beam at 351.1 nm (single line, 100 mW) 
is produced by an argon ion laser (Coherent, Sabre, model DBW25/7), and directed to the first BBO (beta-barium 
borate) crystal (cut for type-I phase matching, optic axis cut at 6 = 35.0°, 1 mm thick) after passing through a 
polarizing beam splitter (PBS) to give a pure horizontal polarization state. Two entangled photon beams (polarized 
in vertical direction) produced by the BBO crystal passing through a sequence of quartz plates, with an optic axis set 
in the diagonal direction. It is pointed out that the photon polarization state will decohere in such an environment 
[ 1 6 1 . The thickness of quartz plates is set in such a way that the photon's polarization decoheres completely (See Ref. 
[27|, the difference between optical path of horizontal and vertical optical polarizations in the quartz plates is 153Ao, 
and Ao is the central wave length of the down-converted light.). The remained pump beam passes through the BBO 
crystal and is rotated to (\H) + \ V)) by a half- wave plate (HWP) and is transmitted through the second BBO 

crystal which is the same as that proposed by Kwiat in Ref. []l7|l8| to produce polarization-entangled photons in 



such a state. The down-converted light beams produced by the first BBO were reflected by four reflectors and passed 
through the second BBO crystal to be mixed with the down-converted beams produced in the second BBO. Note the 
reason that the pump was set in single line instead of single frequency is to decrease the coherent length (about 4 cm 
in our experiment) of the pump light and avoid the interference between the down-converted light beams from the 
two BBO crystals. 

Figure 1. 



At the end, the bi-photon polarization state is measured tomgraphically |L9 18 by quarter- wave plates (QWP), 
HWP and PBS. We use 16 analyzer settings as listed in Table I. The photons are detected by using silicon avalanche 
photodiodes (EG&G, SPCM-AQR) operated in the geiger mode. Each detector is preceded by a small iris (the 
diameter is 1.5 mm) to define the spatial mode, a narrowband interference filter (IF) centered at 702 nm (Andover, 
050FC46-25/7022 [|23f , full width at half is equal to 4.62 nm) to reduce background and define the bandwidth of the 
photons, and a collection lens. The detector outputs are recorded in coincidence with a time-to-amplitude converter 
and a single-channel analyzer, leading to an effective coincidence window of 5 ns. The resulting rate of accidental 
coincidences is less than Is -1 , which can be neglected compared with the typical rate of the true coincidences, which 
is about 300s _1 . 

We also test whether the states produced by this set-up violate the CHSH inequality. This inequality shows that 
15*1 < 2 for any local realistic theory, where 



and E (9\, 9 2 ) is given by 



' = E (6t, 9 2 ) + E (9[, 9 2 ) + E (6 u e' 2 ) - E (9[,0' 2 ) (3) 

c (e 1 , e 2 ) + c (ej ,o£)-c (ej ,e 2 )-c (fli , e£) 



and C (61, 9 2 ) is the coincidence rate of two detectors when the polarization analyzer angels are set in 9\ and 9 2 . In 
this experiment, we selected the settings: 9 X = -22.5°, 6^ = 67.5°, 9[ = 22.5°, 9' 1 ± = 112.5°, 9 2 = 0°, 0f = 90°, 



45°, 0!j- = 135° 
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III. EXPERIMENTAL RESULTS 



We can adjust the intensity proportion (it can be tuned by adjusting the position of the reflectors) between the 
down-converted light from two BBO crystals to obtain Werner states with different coefficients x. Two output states 
have been produced in this experiment, for one of them the CHSH inequality is violated while the other is not. 

The tomographic results are shown in Table I. 

Table I. 



From the data in the third column of Table I we can obtain the density matrix of the first output state directly, 
however, it is not non-negative definite [[l9). We have used the maximum likelihood estimation to construct a 
non-negative definite density matrix 



Pi 



( 0.4169 
0.0203 - 0.0022i 
0.0094 + 0.0237i 
V -0.3476 - 0.0296i 



0.0203 + 0.0022i 

0.0531 
-0.0122 + 0.05271 
-0.0163 - 0.0005i 



0.0094- 0.0237i 
-0.0122 - 0.0527i 

0.0559 
-0.0134 + 0.0191? 



-0.3476 + 0.0296i 
-0.0163 + 0.0005i 
-0.0134-0.0191i 
0.4741 



(4) 



The fit Werner state 



(where x\ = 0.801 ± 0.005) satisfies 



Pi 



Xi 



|$-)($-| + (l-a?i) 



h 



(5) 



F(p 1) p' 1 ) =Tr 2 

= max-F (p-L,p w ) 



11/2 11/2 
Pi PlPl 



(6) 



where pw is an arbitrary two-bit Werner state and the fidelity |24| F (pi, p^) is equal to 0.932 in this experiment. The 
CHSH correlation value [Eq. (3)] |5| = 2.198 ± 0.004 > 2, (the uncertainty of the HWP is A6 ~ 0.2°) which violates 
the CHSH inequality. The theoretical value of |5| of state p[ is equal to 2.266. 

The tomographic results of the second output are shown in the sixth column of Table I. From these data we can 
obtain the density matrix of the second output state 



P2 



0.3949 -0.0217 + 0.0398i -0.0080 + 0.0024i -0.1657 + 0.0150i 

-0.0217 - 0.0398? 0.1285 0.0165 - O.OOlli -0.0075 + 0.0504i 

-0.0080 - 0.0024i 0.0165 + 0.0011* 0.1311 -0.0152 + 0.0228i 

-0.1657- 0.0150i -0.0075 - 0.0504i -0.0152 - 0.0228i 0.3455 



(7) 



The fit Werner state is 



P2 = x 2 |* ) 



(l-*a)£, 



(8) 



where x 2 = 0.405+0.005 and the fidelity F (p 2 , p' 2 ) is equal to 0.982. The CHSH correlation value \S\ = 1.380+0.008 < 
2, which satisfies the CHSH inequality. The theoretical value of \S\ of state p' 2 is equal to 1.146. We should note that 
there is significant difference in count rates between the two sets of results. In our experiment, the intensity proportion 
of down-converted light beams between two BBO crystals is tuned by adjusting the position of the reflectors. So the 
length of down-converted photon's optical path between the source and the detector has been changed for preparing 
the second state. Therefore, the count rates have been decreased significantly. 

We can use the "linear entropy" to quantify the degree of mixture of a quantum state and use the tangle to measure 
the degree of entanglement of a two-bit quantum state |19|,|25|] . The linear entropy for a two-bit system is defined by 



P=l(l~Tr{p*}) 



(9) 



The tangle is defined as T = C 2 , where C is the concurrence in Ref. |25|. In this experiment, for the first density 
matrix in Eq. (4), P = 0.46 + 0.03 and T = 0.35 + 0.01, and for the second density matrix in Eq. (7), P = 0.83 + 0.03 
and T = 0.01 ±0.00. 
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IV. DISCUSSION AND CONCLUSION 



We have prepared two typical Werner states in this experiment and measured the density matrix tomographically. 
Two independent BBO crystals are used to produce parametric down-conversion light beams which can be mixed to 
prepare a wide range of two-bit quantum states. In our experiment, the intensity proportion of down-converted light 
beams between two BBO crystals are tuned by adjusting the position of the reflectors. In fact, the proportion can be 
tuned with linear optical elements, which is shown in Fig. 2. 



Figure 2. 

Although the output light beams were produced by classically mixing (not quantum superposing) the down- 
converted light beams, it is hard (but not in principle) to distinguish from which crystals the output photon is 
produced since the pump light is continuous. However, if we use the pulse pump we can distinguish which crystal 
the output photon from by measuring the photon's arrival time, that is, we can distill the classical information of the 
mixed state and it can be regarded as a pure state. 

We expect that this method can be modified to prepare arbitrary two-bit states conveniently and this source can 
be used in experimental investigation of entanglement purification or other quantum information tasks. 



ACKNOWLEDGMENTS 



This work was supported by the National Fundamental Research Program (2001CB309300), the National Natural 
Science Foundation of China, the Innovation Funds from Chinese Academy of Sciences (CAS), and was also funded 
by the outstanding Ph. D. thesis award and the CAS's talented scientist award rewarded to Lu-Ming Duan. 

Appendix 

For a single photon, the non-dissipative coupling between photon frequency and polarization in a birefringent media 
leads to decoherence of polarization, provided we regard the freedom of frequency as "environment" |l(],[lE],|2(| . 
The single-photon's initial state is described by a pure product state of polarization and frequency, 

|*) = (01 \H)+a 2 \V))® JduiA(u)\u). (A-l) 

with basis \H) (horizontal polarization) and \V) (vertical polarization) denoted by 

|a>-(J).|v>-(! 

respectively, and A (u>) is the complex amplitude corresponding to u>, normalized so that 

doj\A{oj)\ 2 = 1. 

After being transmitted by the birefringent medium whose length is x, the final state of the photon is 

|*' (a)) = ai \H) ® / duA (u) e -fa>(»B-nv)*/c ^ + fl2 ^ g, / du:A ^ e ^(n H -n v ) x /c ^ ( A _ 2 j 



Where njj (nv) is the refractive index of horizontal (vertical) polarization. It is necessary to note that although uh 
and nv are dependent on frequency ui, the value of nji ~~ ny does not vary obviously with frequency varying in a 
small scale and we regard it as constant. The density operator of the polarization state can be obtained by tracing 
over frequency degrees of freedom from the complete density operator 

P(*)=fS5r.^^ r(a) V (A-3) 



ala^T* (x) aia\ 
where 

r (X) = { (LjA (u) A* ( W ) e -Mn H -n v )x/c_ ( A _ 4 ) 
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From Eq. (A-3) it can be seen that the evolution of the polarization's states depend on the character of the 
"environment" which is described by the function A (to). 

The spectrum of the photon's frequency is defined by the narrowband interference filter (IF). In this experiment, 
the spectrum of the photon's frequency is rectangular function. To obtain this environment, we use an interference 
filter whose spectrum's shape is approximate to rectangle and the full width at half maximum (FWHM) is 6 = 4.62 
nm. (Andover, 050FC46-25/7022 @). The theoretical value of T (x) is 

X (TlH — ny) Auj 2c 

and 

x (njj — ny) Auj 



|r(x)| 



sm ■ 



x {n,H — ny) Aui 
where we assume that uh > ny. 

Figure 3. 



2c 



(A-5) 



In the experiment 27 (See Fig. 3, the BBO crystal in the figure is cut at a degenerate type-I phase matching 



angle), we first measure the density matrix of the decohered photon's polarization tomgraphically 

ai a,2 



p = , 

The experimental value of |T (x)\ can be calculated as 

|r(a0| = |a 2 /VaxO3|. (A-6) 

The results are shown in Fig. 4. The coherence of the photon's polarization degenerates as theoretical expected.. It 
is interesting that the coherence is recovered to some extent after it degenerated completely, which is similar to the 
Fraunhofer diffraction at a rectangular aperture (or slit) (26). However, since the spectrum of the filter is not strictly 
rectangular, there are some deviations at the tail of curve. 

Figure 4. 
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Figure caption 

Figure 1. Experimental set-up to produce bi-photon Werner state. 

Figure 2. Experimental set-up to tune the intensity proportion of down-conversion light beams from two BBO 
crystals. 

Figure 3. Experimental setup of controlled decohcrcncc of single photon's polarization. The 351.1 nm line of an 
argon ion laser (Coherent, Sabre) is used to pump a BBO crystal (3 mm thick), which is cut at a degenerate type-I 
phase matching angle to produce a pair of entangled photons. Each photon of the entangled pair is polarized in state 
\H). Then the polarization of one photon of the entangled pair is rotated to (\H) + \V)) by the half- wave plate 
(HWP). After decohering in the quartz plates, this photon's polarization state is measured tomographically by HWP, 
quarter- wave plate (QWP) and polarization beam splitter (PBS). The other photon of the pair is used to do coincidence 
detection. We measure the decohered photon's polarization with 4 analysis settings (\H) , |V) , \H) + \V) , \H) +i \V)), 
allowing reconstruction of the density matrix. 

Figure 4. Experimental results vs theoretical curve. The FWHM of the IF is 4.62nm and the spectrum of the 
frequency is assumed to be rectangular, x is the optical path difference between horizontal and vertical polarizations 
and Ao is the central wavelength — 702. 2nm — of photons produced by SPDC. 

Table I. Settings for measuring bi-photon Stokes parameters. H, V, and D are horizontal, vertical, and diagonal 
(45°) linear polarization, respectively. R is right circular polarization. The data are for state p\ (counted over 100 s). 

Table II. The data are for state pi (counted over 100 s). 
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Fig. 4 



Analyzer Analyzer Coinc. Analyzer Analyzer Coinc. 
in arm 1 in arm 2 count in arm 1 in arm 2 count 



H 


H 


31369 


H 


V 


3945 


H 


D 


18528 


H 


L 


18090 


v 


H 


4180 


V 


V 


35625 


V 


D 


18312 


V 


L 


18557 


D 


H 


17745 


D 


V 


17872 


D 


D 


5257 


D 


L 


21335 


R 


H 


19677 


R 


V 


20088 


R 


D 


20565 


R 


L 


5930 



ti 


u 
ti 


1 n 1 ic* 
iUi Jo 


u 
ti 


V 




ti 


L) 


A 1 


u 
ti 


K 


J05O 


V 


H 


3365 


V 


V 


8867 


V 


D 


5725 


V 


R 


5532 


D 


H 


6546 


D 


V 


5891 


D 


D 


3829 


D 


R 


5210 


R 


H 


6688 


R 


V 


4788 


R 


D 


5085 


R 


R 


7272 



Table I 



